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ABSTRACT 


This paper presents an approach for acquiring a magnetic particle imaging 
(MPI) signal, by utilizing the second harmonic detection of the magnetic 
nanoparticles tracers. An MPI signal with high signal-to-noise ratio (SNR) is 
crucial for high spatial resolution images that reveals the distribution of the 
tracers in a target area. Samples of Resovist and Perimag nanoparticles 


tracers were prepared in liquid and immobilised form, which were placed at 

some distances under the receiver coil of the single-sided MPI scanner. 
Keywords: The samples were exposed to the excitation magnetic field generated at 22.8 
kHz and a static gradient field generated with a direct current of 2 A. 
The non-linear magnetization response of the tracers for each spatial position 
is recorded in the form of voltage signal by a gradiometer pickup coil, 
with the second harmonic signal being extracted by a resonance circuit. 
The results obtained revealed that a sufficient signal from the tracers is 
recorded at up to 25 mm under the pickup coil, with Perimag samples 
inducing higher signals as compared to Resovist. The dependence of the DC 
gradient field on the second harmonic signal shows that the peak signal 
amplitude for Resovist and Perimag particles as +5 mT and +6 mT 
respectively. Additionally, the second harmonic signal amplitude increases 
exponentially with an increase in the excitation magnetic field. 
Thus, the results obtained shows the potential of this approach in acquiring 
high SNR MPI signals at low excitation frequency, which could be vital in 
reconstructing the contour images of the tracers, particularly in sentinel 
lymph node biopsy (SLNB) for breast cancer diagnosis. 


Magnetic nanoparticles 
Magnetic particle imaging 
Second harmonic detection 
Tomographic images 


Copyright © 2019 Institute of Advanced Engineering and Science. 
All rights reserved. 


Corresponding Author: 


Muhammad Mahadi Abdul Jamil, 
Department of Electronics Engineering, 
Universiti Tun Hussein Onn Malaysia, 
Johor, Malaysia. 

Email: mahadi@uthm.edu.my 


1, INTRODUCTION 

Magnetic particle imaging (MPI) is a new imaging method which exploits the non-linear 
magnetization response of magnetic nanoparticles (MNPs) tracers to reveal their distribution within a 
target [1]. The MNPs are excited by the AC magnetic modulation field generated by the MPI scanner and 
respond with a non-linear magnetization. The magnetization (usually called MPI signal), contains not only 
the fundamental excitation frequency, but also its harmonics occurring at multiples of the excitation 
frequencies [2-5]. Harmonic detection of MNPs is a non-linear AC susceptibility detection approach, where a 
sufficient excitation field with amplitude HO is applied to magnetize the MNPs. The influence of the 
excitation field 1s then removed by a Bandstop filter, and the amplitude of the harmonics from the non-linear 
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magnetization response of the particles is utilized for the MNPs image reconstruction. Before image 
reconstruction, the fast Fourier transform (FFT) of the signal is determined, to convert the signal from time 
domain to frequency domain, where the second harmonic signal is certainly picked by a resonant circuit and 
logged into the computer by a data acquisition device (DAQ) [3, 6]. 

In this paper, an alternative approach of acquiring the second harmonic MPI signal is proposed and 
demonstrated experimentally. Initially, a gradiometer coil sensor is employed as the MPI signal receiver in 
order to cancel out the interference of the excitation field and the earth’s magnetic field. Then, a Bandstop 
filter is used to further attenuate the component of the excitation, so as to have a sufficient signal with high 
signal-to-noise ratio (SNR). Finally, the second harmonic signal amplitudes are extracted by the resonant 
circuit designed for second harmonic detection. These steps of MPI signal acquisition yields a high SNR 
signal which increases the spatial resolution of the reconstructed images, which could help in sentinel lymph 
node biopsy (SLNB) for breast cancer diagnosis. 


2, DETECTION METHODS OF MNPS 
Thus far, there are three techniques for MNPs detection as reported in [7], which are: AC 
susceptibility (ACS), magnetic relaxation (MRX) and remanence-based measurement methods. In ACS, a 


periodic excitation field y =H 9 sin(2af yt) having an excitation frequency fo and amplitude H 0 1s applied 


to the MNPs. When the value of Hy is small such that the Langevin parameter is less than unity, 
the magnetisation produced by the MNPs is within the linear region. This situation is called linear ACS, and 


the signal frequency cannot be differentiated from the excitation frequency. But when the Hp 1s large, 


a weak signal is obtained, that needs to be extracted from the large Hy). The extraction process utilizes the 


non-linear magnetisation response of the MNPs. This condition is called non-linear ACS, where the 
harmonics generated due to the non-linear magnetisation of the MNPs are employed, and easily differentiated 
from the excitation frequency. There are two types of non-linear ACS; harmonic detection approach and the 
mixed-frequency detection approach [8-10]. In harmonic detection approach, a sufficient is used to fulfil the 
condition that the Langevin parameter is above two (¢ > 2), in which even the third harmonic can be 
identified [11]. Yoshida et al., (2012) developed a system for the detection of magnetic nanoparticles using 
the second-harmonic signal from the MNPs, using a simple inductive pickup coil gradiometer [12]. 
The advantage of using the second-harmonic signal is that the interference of the excitation field is 
significantly reduced. Furthermore, [13] developed an MPI system that uses harmonic signals of the MNPs 


by their non-linear magnetisation curve. Harmonic signals from the MNPs were measured when wpH gq, of | 
mT is applied to the sample while “9H ,. was varied from 0 to 5 mT, and a solid-phase MNPs ware used. 
Mixed-frequency detection approach, on the other hand, uses two different excitation fields Ho, 


and Hy, , with H ¢, having low frequency but with high amplitude and H_, having high frequency but with 


small amplitude of excitation. Therefore, a response containing two frequencies 1s produced by the MNPs. 
This response is used in identifying the MNPs [14]. In MRX, the MNPs are exposed to a high magnetic 
excitation after an initial time ¢. H is switched off and the decay in the magnetic signal generated by the 
MNPs (relaxation signal) is measured for +>0 using another magnetic field (pulse magnetic field) [15]. 
Magnetic relaxometry was used to detect breast cancer cells, using targeted MNPs and ultra-sensitive 
magnetic field sensors, in particular using superconducting quantum interference device (SQUID) sensors 
in [16]. The magnetic remanence-based approach is comparable with MRX, as MNPs with long Neel 
relaxation period are employed. Here, the MNPs response is retained for an elongated period to differentiate 
the positions at which the excitation and the signal measurement are made [17]. In this study, the harmonic 
detection approach is implemented in order to differentiate the harmonics easily by applying a sufficient 
excitation field amplitude. 


3. RESEARCH METHODS 

For the second harmonic signal detection to be conducted, an MPI scanner is required. For that 
reason, a single-sided MPI scanner developed by [18] is used in this study. The scanner is a one-dimensional 
scanner having 22.8 kHz excitation frequency. In addition, Resovist (Fuji pharma) and Perimag (Micromod) 
MNPs were employed as the tracer materials, both immobilised and in liquid form to imitate 
their accumulation in the tumour region in real life SLNB. Explaining research chronological, including 
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research design, research procedure (in the form of algorithms, Pseudocode or other), how to test and data 
acquisition [1-3]. The description of the course of research should be supported references, so the explanation 
can be accepted scientifically [2, 4]. 


3.1. Single-Sided MPI Scanner Set-Up 

The measurement system for the second harmonic detection of MNPs is depicted in Figure |, which 
consists of the constructed coils, arranged coaxially on the same plane to form the single-sided MPI scanner. 
The designed BPF and BSF, as well as the resonant circuit (for second harmonic detection) were also 
connected appropriately to the transmitting and receiving sections of the scanner device [19, 20, 21, 22]. 
Other equipment utilised in the scanner realisation is also shown in Figure 1. These equipment are the 
function generator (BK Precision 4055), AC power amplifier (HP 8447F H64), DC power supply (Tektronix 
PS280) and the DAQ device (NI USB-6009). Additionally, sheets of acrylic material (9 mm thickness) were 
also used to hold the coils coaxial arrangement in a fixed position. The acrylic material is selected due to its 
little weight and does not exert any magnetic field on the scanner operation. Furthermore, it has a reasonable 
degree of compatibility with human tissues. 

















Figure 1. Measurement system set-up 


The excitation field is generated by an alternating current of 1 A, with an excitation frequency of 
22.8 kHz, having an amplitude of 30 mV peak-to-peak. Conversely, the selection field is generated with a 
direct current of 2 A. The currents for the excitation and selection field were applied in the opposite direction, 
to achieve spatial encoding and the creation of the field free point. The parameters of the coils that made up 
the single-sided MPI scanner are summarised in Table 1. 


Table 1. Coils Specifications 


Tvpe Excitation Selection Pickup 
YP Coil Field Coil Coil 
Number of turns 352 416 208+208 
Number of layers 8 9 8 
Wire diameter 06 06 0.6 


(mm) 


3.2. Single-Sided MPI Scanner Set-Up 

As explained earlier, two commercially available MNPs were employed as the tracer materials, to be 
detected under the influence of the magnetic fields generated by the developed MPI scanner. Therefore, 
samples of Resovist (Fuji pharma) and Perimag (Micromod) particles were prepared for the second harmonic 
detection and inserted into one well of a 96-well cell culture plate (655-180, flat-bottom), with a total sample 
volume of 300 ul assumed. The iron concentration in the MNPs found at the sentinel lymph node is estimated 
to be around 100 ug [2]. Thus, 100 ug Fe is employed as the iron concentration in the MNPs samples to be 
detected using second harmonic detection. 

Initially, the samples were diluted with distilled water to achieve the desired volume in liquid form 
and mixed with glycerol to get the immobilized phase. Consequently, the volume of the MNPs needed to 
contain 100 ug Fe is calculated using the particles parameters given in Table 2, as 3.6 ul for Resovist and 4.0 
ul for Perimag. Finally, a volume of 300 ul sample is attained by adding 296.4 ul of distilled water to 
Resovist and 296 ul of distilled water to Perimag MNPs. 
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Table 2. Resovist and Perimag MNPs Parameters 


MNPs Concentration Density Calculated volume 
(mg Fe/ml) (g/ml) (ul) 
Resovist 
(60nm) 28 1.057 3.6 
Perimag 
(130nm) 25 1.400 4.0 


The MNPs need to be experimented not only in fluid form but also in a solid phase ((¢mmobilised). 
The samples are required to be immobilised to imitate the interaction between the tissues and the MNPs as in 
real life scenario. In this regard, the MNPs suspensions were diluted with glycerol to prepare the immobilised 
samples. The glycerol serves as a coagulant, which changes the fluid into semi-solid form. Hence, Resovist 
and Perimag MNPs were mixed with 235 ul and 235.2 ul of glycerol respectively, to form an immobilised 
sample as shown in Figure 2. 





Figure 2. Immobilised MNPs samples (a) one sample perimag and (b) two sample resovist particles 


3.3. Experiment Performed 

The experiments to detect the tracer samples in liquid and solid form using second harmonic 
detection were conducted. For proper investigation on the harmonious features of the MNPs samples, 
the excitation magnetic field was fixed at 0.8, 1.7, 2.6, 3.0, 3.4 and 4.4 mT at 22.8 kHz excitation frequency. 
Besides, the selection magnetic field was varied from 0 to 5 mT. Initially, the samples were placed under the 
scanner (one at a time), exposing them to the excitation magnetic field at the FFP by moving the sample to all 
the spatial positions across the FOV. Since the developed MPI scanner is one-dimensional, the FOV has a 
dimension of 60 x 60 mm? with a total of 169 spatial locations. Due to the influence of the excitation field, 
the MNPs samples became saturated and responded by producing another magnetic field. The particle’s 
magnetisation response for each spatial position is recorded in the form of the voltage signal, at the second 
harmonic frequency. Measurements of the particles magnetisation were executed in the z-direction, where the 
excitation magnetic field is assumed to be highly uniform in strength and path. The right-hand rule for 
magnetic field produced by a solenoid states that the thumb points in the direction of the field, while the 
coiled fingers indicate the direction of current [23]. The magnetic field diminishes with increasing distance 
away from the coil. Therefore, the distance along the z-direction is varied to determine how far the second 
harmonic signals can be detected, which is related to the distance between the MPI scanner and the sentinel 
lymph node under examination. In this regard, the mean depth of the sentinel lymph node from the skin 
surface in humans is estimated to be 12 + 5 mm [2, 24]. Hence, the distance in the z-direction is varied from 
the surface of the pickup coil to about 30 mm, so that the sentinel lymph node distance is appropriately 
covered. Lastly, the voltage induced in the pickup coil at each interval in the z-direction is recorded, with the 
second harmonic signal amplitude determined through FFT of the induced voltage signal, which is employed 
in the MNPs image reconstruction. 


4. RESULTS AND ANALYSIS 

Initially, the MNPs samples were placed (one after the other) at the centre of the designed single- 
sided MPI scanner for excitation and the corresponding MPI signal measurement. The vertical 
distance (z-direction) between the sample and the scanner is varied from the surface of the pickup coil to over 
30 mm in space. The variation in the z-direction determines the sensitivity of the developed scanner with 
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increasing distance. Liquid samples of Resovist and Perimag particles were prepared, each containing 100 ug 
of the tracer material in a 300 ul volume. Immobilised samples, on the other hand, contains 235 ul of 
glycerol and the same tracer content as the liquid sample (100 ug). The excitation magnetic field was 
generated with a periodic signal of 22.8 kHz, from an AC of | A while the DC for the selection field was set 
at 2 A. As a result, the induced voltage at the second harmonic when the distance in the z-direction is varied 
is depicted in Figure 3. 

Perimag particles induced higher voltage signals as compared to Resovist particles up to a distance 
of 25 mm, as shown in Figure 3. Furthermore, liquid samples produce greater voltage signals than the 
immobilised samples, which is due to more rapid decay of Neel relaxation times and the size difference of 
the MNPs as explained earlier. Thus, both Resovist and Perimag particles could be evidently detected at 25 
mm in the z-direction, which is adequate for tracer detection in SLNB application. 


4.1. Second Harmonic Signal Measurements 

Subsequently, the MNPs tracer sample is scanned through the FOV and the amplitude of the second 
harmonic induced signal is recorded, to determine the position of the MNPs sample in the XY-plane by 
plotting its contour map. Before the contour plot, the dependence of the second harmonic signal with varying 
the DC magnetic field for the selection field (Hac) 1s investigated. The Hac 1s varied between +25 mT to have 
signals induced from both the positive and negative magnetisation of the MNPs, while the AC excitation 
field is fixed at 10 mV peak-to-peak from 22.8 kHz. Figure 4 depicts the second harmonic signals recorded 
from varying Hdc for Resovist sample. It can be deduced that the signal amplitude for liquid Resovist sample 
reaches its peak of 15 mV at around 5 mT Hdc, while the immobilised sample has a peak signal amplitude of 
approximately 13 mV at the same DC magnetic field. 





—e— Resovist-Solid 
—e— Resovist-Liquid 
—— Perimag-Solid 
—— Perimag-Liquid 


2.5x10% 4 | 











2.0x10% 4 + 








1.5x10* 


1.0x10* 


Induced voltage (V) 


5.0x10° 





0.0 
) 10 20 30 40 
Distance, z (mm) 


Figure 3. Induced voltage with varying distance along z-direction 
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Figure 4. Second harmonic signal dependence on different DC magnetic field for Resovist sample 


The second harmonic signal dependence with varying DC magnetic field for Perimag sample was 
also investigated and is presented in Figure 5. Liquid Perimag sample generates peak signal amplitude of 
around 15.8 mV as compared to 14 mV produced by the immobilised Perimag particles at +6 mT DC 
magnetic field. For symmetry and compliance with the AC excitation magnetic field applied in both negative 
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and positive sides of the particles magnetisation curve, the DC magnetic field 1s also applied in a negative 
and positive direction, to have signals induced in both directions. 

It is therefore evident from the measurements conducted on the four samples that the dynamic 
region of Resovist and Perimag particle magnetisation curve is within +5 mT DC magnetic field for Resovist 
and 6 mT DC magnetic field for Perimag, after which no further magnetisation is recorded, as revealed by 
Figure 4 and Figure 5. This implies that the magnetisation response of the tracers is confined to the FFP by 
the selection field at +5 mT for Resovist and +6 mT for Perimag particles. Additionally, the voltage signal 
induced by the tracers is purely contributed by the tracers within these DC magnetic fields. 

Next, the dependence of the AC excitation magnetic field on the induced signal is explored, by 
varying the excitation field between O-10 mT while the DC selection field was kept constant at 5 mT. 
The results obtained revealed that the signal amplitude of the second harmonic frequency increases 
exponentially with increasing Hac, as shown in Figure 6. For Resovist particles, similar signal amplitudes 
were recorded for Hac < 2 mT, as a result of similar averaged relaxation time experienced by the liquid and 


immobilised samples, in addition to being under the same temperature condition. 
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Figure 5. Second harmonic signal dependence on different DC magnetic field for Perimag sample 
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Figure 6. Dependence of excitation magnetic field on the induced signal for resovist 


Perimag particles, on the other hand, records different signal amplitude even at low excitation field, 
with the liquid sample generating higher harmonics than the immobilised sample as depicted in Figure 7. 
The various signal amplitudes obtained in the measurement is as a result of faster Neel relaxation time of the 
liquid sample as compared to that of the immobilised sample, due to the mean time between two flips of the 
particle’s internal magnetisation. The signal amplitudes recorded for Perimag samples show a similar trend to 
that obtained from Resovist samples, with the exponential rise [25, 26]. 

The dependence of DC magnetic field on the excitation magnetic field for the peak second 
harmonic signal detection is finally investigated. As can be seen in Figure 4 and Figure 5, zero second 
harmonic signal amplitude was recorded when Hac is zero for both Resovist and Perimag particles but reaches 
its peak value when Hac increases. Figure 8 revealed a virtually linear relationship between Hac and Hac, as 
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the small value of Hac peak could be obtained with little Ha. when the highest amount of the signal amplitude 
is recorded [12, 13, 26-28]. 
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Figure 7. Dependence of second harmonic signal amplitude on AC excitation magnetic field for Perimag 
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Figure 8. Relationship between Hac and Hac at peak signal amplitude 


5. CONCLUSION 

It has been shown that MPI signal with high SNR could be appropriately acquired from the non- 
linear magnetization response of the tracers, using low frequency excitation magnetic field. The use of 
Bandstop filter removes the component of the excitation field, yielding a high SNR signal completely from 
tracer materials. Additionally, the simple resonant circuit designed for the second harmonic signal detection 
was found to be appropriate in extracting the MPI signal at the second harmonic frequency. Measurement 
results obtained clearly revealed the width of the dynamic range of Resovist and Perimag nanoparticles, with 
the latter producing higher signal amplitudes in all experiments, due to its large hydrodynamic diameter. 
Liquid samples were found to produce higher signal amplitudes than the immobilised samples due the 
dominance of Brownian rotation over Neel rotation of the particles. Thus, the outcome of this research 
showed that the MPI signal acquisition approach has a potential in the detection of tracers, which could help 
in obtaining high resolution tomographic images, in sentinel lymph node biopsy (SLNB) for breast cancer 
diagnosis. 
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